Prioritizr Runs

We used the fPrioritizrRun() function to streamline the runs. fPrioritizrRun() arranges all the data and makes sure that they’re in the right format, then runs problem() and solve(). This function also saves the summaries each plan (cost, % area protected and number of representative features that achieved 30% or more protection).

The relevant scripts for this step is found in ~/scripts/ with the prefix 10_Prioritizr.

# the function is in:
source("scrpits/10a_Prioritizr_fPrioritizr.R")
# to run the function with climate-smart data
source("scripts/10h_Prioritizr_fPrioritizrIUCNSmart_NoProv.R")
# to run the function with climate-uninformed data
source("scripts/10i_Prioritizr_fPrioritizrIUCNUninformed_NoProv.R")
# to create no-regret closures with the climate-smart plans

To create no-regret plans, we used thre function fCreateNoRegret() which intersected the selected PUs in all three climate scenarios. The relevant scripts for this step is found in ~/scripts with the prefix 11_NoRegret.

# the function is in:
source("scripts/11_NoRegret_fCreateNoRegret.R")
# the runs are found in:
source("scripts/11e_NoRegret_fCreateNoRegretIUCNRuns_NoProv.R")

As mentioned in the previous section, we will only be showing results of Runs 2 (Max = 0.9, Min = 0.1), 5 (Max = 0.6, Min = 0.1) and 8 (Max = 0.3, Min = 0.1).

Run 2: Maximum Target = 90%

Climate-smart runs

SSP1-2.6 run

SSP126_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP2-4.5 run

SSP245_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP5-8.5 run

SSP585_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

No-regret plan

noregret_target90 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target90.rds")

Climate-uninformed run

uninformed_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 12
##   target_scenarios       sum_area percent_area num_pu total_cost median_velocity
##   <chr>                     <dbl>        <dbl>  <dbl>      <dbl>           <dbl>
## 1 Target90_summaryuninf…   6.72e7         79.0  25207    767309.       67242193.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## #   median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## #   scenario <chr>

Run 5: Maximum Target = 60%

Climate-smart runs

SSP1-2.6 run

SSP126_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP2-4.5 run

SSP245_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP5-8.5 run

SSP585_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

No-regret plan

noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")

Climate-uninformed run

uninformed_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 12
##   target_scenarios       sum_area percent_area num_pu total_cost median_velocity
##   <chr>                     <dbl>        <dbl>  <dbl>      <dbl>           <dbl>
## 1 Target60_summaryuninf…   4.48e7         52.7  16807    393372.       44834353.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## #   median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## #   scenario <chr>

Run 8: Maximum Target = 30%

Climate-smart runs

SSP1-2.6 run

SSP126_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP2-4.5 run

SSP245_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

SSP5-8.5 run

SSP585_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## #   percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## #   median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## #   represented_features <dbl>, target <dbl>, scenario <chr>

No-regret plan

noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")

Climate-uninformed run

uninformed_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 12
##   target_scenarios       sum_area percent_area num_pu total_cost median_velocity
##   <chr>                     <dbl>        <dbl>  <dbl>      <dbl>           <dbl>
## 1 Target30_summaryuninf…   2.24e7         26.3   8407    149748.       22426513.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## #   median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## #   scenario <chr>